Gaussian wave packets in phase space: The Fermi gp function 
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Any pure quantum state can be equivalently represented by means of its wave function ip{q) or 
of the Fermi function gF{q,p), with q and p coordinates and conjugate momenta of the system 
under investigation. We show that a Gaussian wave packet can be conveniently visualized in phase 
space by means of the curve = 0. The evolution in time of the gp — curve is then 

computed for a Gaussian packet evolving freely or under a constant or a harmonic force. As a 
result, the spreading or shrinking of the packet is easily interpreted in phase space. Finally, we 
discuss a gedanken prism microscope experiment for measuring the position-momentum correlation. 
This gedanken experiment, together with the well-known Heisenberg microscope and von Neumann 
velocimeter, is sufficient to fully determine the state of a Gaussian packet. 



I. INTRODUCTION 

The free propagation of a Gaussian wave packet is the 
simplest application of the time-dependent Schrodinger 
equation. In their first quantum mechanics course, stu- 
dents soon learn that, for a free particle, a minimum- 
uncertainty Gaussian packet spreads in time. They also 
learn that, for a harmonic oscillator, such minimum 
packet does not spread. Since, according to the Ehren- 
fest theorem, a quantum wave packet follows a beam of 
classical orbits as long as the packet remains narrow, one 
can conclude that the harmonic motion of a minimum 
packet can be described by means of classical mechanics. 

On the other hand, it should be stressed that: 

(i) If we call q and p the coordinate and momentum 
of the system under investigation, then also a beam of 
classical orbits can spread along q. That is, the spreading 
is not necessarily a quantum phenomenon. 

(ii) Gaussian wave packets can also contract along 
g, namely the position uncertainty can narrow along q. 
Given an arbitrarily long time t, it is possible to build 
packets that shrink up to that timc.^ 

Heisenberg uncertainty principle forbids the notion of 
trajectories in quantum mechanics. However, a com- 
parison between classical and quantum dynamics can be 
made by studying the evolution in time of the classical 
and quantum phase space distributions. The 
classical phase space distribution pc {q, p, t) determines 
the probability pc{q,P,t) dqdp that a particle is found 
at time t in the infinitesimal phase space volume dq dp 
centered at {q,p). The evolution in time of p is governed 
by the Liouville's equation.^ A convenient phase space 
formulation of quantum mechanics is obtained using the 
Wigner distribution function p\y{q,P,t).^'^'^^ Such distri- 
bution evolves according to a quantum Liouville equa- 
tion. However, the Wigner function in general is not 
positive definite and therefore requires the introduction 
of delicate concepts such as that of quasiprobability. 

In this paper, we follow a different phase space ap- 
proach, based on an old paper by Fermi. The relevant 



quantity here is the Fermi gp function, or more precisely 
the gp{q,p,t) = curve. We will show that the gp func- 
tion provides a nice intuitive phase space picture of Gaus- 
sian wave packets and of their evolution in time. In par- 
ticular, the spreading or shrinking of a Gaussian packet 
evolving freely or under a constant or a harmonic force 
are easily interpreted in the Fermi-function picture. Fi- 
nally, we will discuss gedanken experiments determining 
the gp ~ Q Fermi curve for Gaussian packets. In particu- 
lar, we will focus on a prism microscope capable of mea- 
suring the position-momentum correlation. It turns out 
that the prism microscope, together with the well-known 
Heisenberg microscope and von Neumann's velocimeter, 
fully determines the gp = Q curve for Gaussian packets. 

II. THE FERMI gp FUNCTION 

For the sake of simplicity, we consider the case of a 
single particle moving along a straight line. In classical 
mechanics, it is possible to determine the state of such 
system by measuring at some time t the values of the 
position q and momentum p or, equivalently, by measur- 
ing two independent functions f{q,p) and g{q,p), from 
which the values of q and p can be obtained. 

In quantum mechanics, due to Heisenberg uncertainty 
principle, a simultaneous exact measurement of q and p is 
not possible. An exact measurement is in principle possi- 
ble only for one physical quantity: q or p or any function 
g{q,p)- Of course different choices oi g{q,p) lead to differ- 
ent determinations of the system. For instance, an exact 
measurement of q implies a complete indetermination of 
p, and vice versa. 

As pointed out by Fermi, the state of a quantum 
system at a given time t may be defined in two com- 
pletely equivalent ways: (i) by its wave function i})[q,t) 
or (ii) by measuring a physical quantity gp{q,p,t). In- 
deed, given the measurement outcome gp{q,p,t) = g, 
then ip{q., t) is obtained as solution of the eigenvalue equa- 
tion gp{q,p,t)ip{<i,t) = g'ip{q,t), where p = —ihdq. On 



the other hand, given the wave function tp{q,t) it is al- 
ways possible to find an operator gF{Q,P,t) such that 



9F{q,P,t)tp{q,t) = 0. 
Using the polar decomposition 

V;(g,t) = p(g,t)e^^(«'*\ 



(1) 



(2) 



where p and 6 are real (p > 0), it is easy to check that 
identity (1) is indeed fulfilled by taking 

gp {q, -ihd„t) = [-ihdg - fid,e{q,t)f + h^^ll^^. 

Z"^' (3) 

Equation (1) implies that the corresponding physical 
quantity gF{q,P,t) takes the value g = 0- Therefore, the 
equation 



gF{q,P,t) = \p- hdg0{q,t)f + h- 



(4) 



defines a curve in the two-dimensional phase space, para- 
mctrically dependent on time t. In other words, as ex- 
pected from Heisenberg uncertainty principle, we cannot 
identify a quantum particle by means of a point {q,p) 
but we need a curve gF{q,p,t) = 0. We call gp the 
Fermi function and we are interested in the phase space 
curve gp = 0. It is also possible to write equation (4) in 
the form 



p± = fu)q9{q, t) ± 



(5) 



This equation locates two points {q,p+) and {q,p-) in 
the phase space for any q such that dgqp{q,t) < and 

P{q,t)^0.'' 



III. THE FERMI qf FUNCTION FOR 
GAUSSIAN WAVE PACKETS 

In this section, we discuss the phase space curve 
gpiq^P, t) = corresponding to the wave function ijj{q,t) 
solution to the one-dimensional Schrodinger equation 



ihdti){q,t) 



2m 



V{q) 



(6) 



for various time- independent potentials V{q) such that 
an initially Gaussian wave packet tl>{q,0) remains Gaus- 
sian at all times. 



A. Free particle 

We have V{q) =0 and consider as initial condition the 
Gaussian minimum uncertainty wave packet 



V(g,0) = 
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FIG. 1: Phase space curves = at different times (from left 
to right, T = = —3, —2, —1, 0, 1, 2, 3) for a free particle, 
whose state at t = is a minimum uncertainty Gaussian wave 
packet centered at = 0, ^ = 2). 



The solution to the Schrodinger equation (6) with this 
initial condition reads 



ilj{q,t) 
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We now apply operator (3) to this wave function, thus 
obtaining 



fit)gF{q,p,t) = ^ + ^ ( 1 



p 



m 



m(52 
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where for later convenience we have multiplied the Fermi 
function by f{t) = ^ (l + j and we have defined 



Pot 

q = q-qo , p = p-po- 

m 

The gF{q,P,t) = curve is the ellipse 
aq^ + bp^ + 2cqp = 1, 

where 
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As shown in Fig. 1, the shape of the ellipse changes over 
time. However, its area 
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remains constant. 



B. Uniformly accelerating particle 

We have V{q) = —F^q. In this case, the solution to 
the Schrodinger equation (6) with initial condition the 
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FIG. 2: Same as in Fig. 1, but for a uniformly accelerating 

particle with = 1.5. The t = Q minimum uncertainty 

Gaussian wave packet is centered at = 0, ^ = O). 



Gaussian minimum uncertainty wave packet (7) reads 

(14) 




where 



' = F^qt - 
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(15) 



The Fermi gp function is again given by equation (9), 
but with the center of the ellipse accelerating uniformly, 
namely we have 



Pot 

qo 

m 



2m ' 



p = p-po- Fot. (16) 



The deformation with time of the ellipse gp = is inde- 
pendent of the force Fq, which simply displaces the center 
{q, p) of the ellipse along the parabola 



9 = ?o 



""-(p-^>o) + ^^. (17) 



m_Fo 2mFQ 
The curves gp = at different times are shown in Fig. 2. 

C. Harmonic oscillator 

Wc have V{q) = ^mtj^q^. A Gaussian solution to the 
Schrodinger equation (6) is given by 

V(g,t) = N{t)exp^^-^A{t)[q-Q{t)]''+ix{q,t)^ , 

(18) 



N{t) = 




C{t) + iVBS{t) 



(19) 



C{t) = cos(a;oi + (f)), S{t) = sm{u;ot + </>), B 

A{t)=Ar.{t)+iAi{t), 
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Ar{t) = 



h[c^{t) + Bs^t)y 
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(20) 
(21) 

(22) 

(23) 

(24) 
(25) 



Different values of the parameters a, Qo and (p corre- 
spond to different solutions to the Schrodinger equation. 
Note that, for the sake of simplicity, we have not consid- 
ered in (18) the most general Gaussian wave packet for 
the harmonic oscillator. 

Due to the time dependence of A(t), it is not immediate 
to grasp the evolution of the Gaussian packet (18). A 
clear illustration is provided by the Fermi function. The 
gp{q,p,t) = curve is derived from equation (4) and 
again is a (paramctrically dependent on t) ellipse in phase 
space, whose equation can be written in the form (11), 
with 



Al + A^, 



b = 



h'^Ar' 



c = 



Q = <l-Q{t), p = p-mQ{t). 



(27) 



Note that the area of the ellipse remains constant: A = 
nh. Owing to the term proportional to qp, the ellipse in 
general does not have its axes parallel to the coordinate 
axes q,p. This term disappears when Ai ~ 0, that is, 
for B = 1 (a^ = H^). This special case corresponds to 
coherent states, for which the gp = curve becomes, in 
the coordinate plane ^q, i^^^ , a rigidly moving circle: 



lp-mQ{t)Y^ 



rrwoo 



(28) 



A more generic example (squeezed state) is shown in 
Fig. 3. 
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case of a free or uniformly accelerating particle we obtain 



q/Qo 



FIG. 3: Same as in Fig. 1, but for the squeezed state (18) of 



a liarmonic oscillator, with 



20, B = 0.1 and woi + < 



ranging from to 77r/4, in steps of 7r/4. 



IV. MEASURING THE 
POSITION-MOMENTUM CORRELATION: A 
GEDANKEN EXPERIMENT 

It is easy to check that, for the examples discussed 
above, the evolution of the phase space region enclosed 
by the gF{Q,P,t) = curve is the same as for the cor- 
responding classical phase space distribution pc{q,p,t), 
uniform inside such region and zero-valued outside and 
evolving according to the Liouville equation. Due to this 
equivalence, the conservation of the area of this region is 
assured by the Liouville theorem.^ However, we should 
stress that the spreading or the shrinking of a wave packet 
are quantum effects, in the sense that the wave packet 
refers to a single particle, while the distribution function 
Pc{Q,P,t) describes a classical ensemble of orbits. No 
physical meaning should be attached to each point in- 
side the gp — curve. Due to Heisenberg principle it is 
not possible to assign with arbitrarily small uncertainties 
both position and momentum. 

For Gaussian states, the area enclosed by the Fermi gp 
function has a simple interpretation in terms of Heisen- 
berg uncertainty principle. To illustrate this point, let us 
compute the variances 

(Aqf = {{q {q})'), {Ap)' = {{p - (p))^) (29) 

and the position-momentum correlation 

K = l{{q- {q)){p - (p)) + {P- mi - m (30) 
for the examples discussed in the previous section. In the 
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while for the harmonic oscillator 
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In both cases, the generalized uncertainty relation^^'-'^'''-'^* 



iAqr{Apr-K' = ^ 



(33) 



is fulfilled. Moreover, the parameters determining the 

gp = ellipse can be expressed as follows: a = ■p-(Ap)^, 
b = -^{Aq)^ and c = -^K (for the free or uniformly 
accelerating particle) or c = (for the harmonic 

oscillator). This implies that in these examples the ex- 
perimental determination of the gp = ellipse is possible 
once the quantities {q), {Aq)^, (p), (Ap)"^ and K are mea- 
sured. 

A simple consequence is that the spreading (in q) of 
a Gaussian wave packet for a free or uniformly acceler- 
ating particle as well as the oscillations in time of the q 
and p variances for a harmonic oscillator can be simply 
visualized in terms of the evolution in time of the gp = 
curve or, equivalently, of the phase space region enclosed 
by this curve. We stress that here the conservation of 
the area of the gp = ellipse is equivalent to the validity 
of the generalized uncertainty relation (33) at all times. 
The product AqAp may grow or even reduce in time. 
This can be clearly seen from Figs. 1 and 2: AqAp grows 
if we start from a minimum uncertainty state; on the 
other hand, if we start from a different state, AqAp may 
decrease until the minimum uncertainty state is reached 
and then starts growing. 

Notwithstanding the behaviour of AqAp, the uncer- 
tainty in the measurement of the state in phase space 
is constant, provided that we measure not only position 
and momentum but also their correlation K. For the 
complete determination of the gp = ellipse we need 
three thought experiments: 

(i) The measurement of (q) and (Aq)'^. This can be 
done by means of Heisenberg microscope. 

(ii) The mesurement of (p) and (Ap)^, possible using 
von Neumann velocimeter. 

(iii) The measurement of the position-momentum cor- 
relator K. 

The first two gedanken experiments are well known and 
their description can be found, for instance, in Ref.^®. In 
order to measure the correlator K, one can use Heisen- 
berg microscope with the addition of a prism after the 
objective, ^° see Fig. 4. 

The main idea of this prism microscope gedanken 
experiment is suggested by the objective prism spec- 
troscopy used for astronomical spectroscopy. In that 



5 



s 




FIG. 4: Schematic drawing of a thought microscope, mea- 
suring the position-momentum correlation of a particle. The 
scattered light is collected on the screen S. Note that a coin- 
cidence system with the detector D of the scattered particle 
is necessary, in order to fix the scattering angle. 

case, a prism is placed in front of the objective lens of 
a telescope. As a result, the light coming from stars 
within the field of view of the telescope is dispersed by 
the prism depending on its frequency. Therefore, the 
telescope can be used for simultaneously measuring the 
spectra of many stars. In our case, the prism is located 
after the lens of the microscope (see Fig. 4), so that, if 
the refractive index of the prism is a suitable function of 
the wave length of the scattered photon, then such pho- 
ton strikes the screen S" in a position dependent of the 
position-momentum correlation of the measured particle. 

We now sketch the working of the prism microscope. 
We first write the equations of the Compton effect for 
the scattering between the photon (our probe) and the 
measured particle. Considering a scattering event on 
the plane of Fig. 4, the nonrelativistic energy- momentum 
conservation laws read as follows: 

{hiyQ + -mc^Pl — hv + 
hvQ — m(3oC^ — hv cos 6 -\- mPc^ cos ip, ^'^^'^ 
hv sin 6 — m(3c^ sin Lp = 0, 

where h is Planck constant, c the speed of light in vac- 
uum, m the particle mass, uq and u (/Jqc and /3c) denote 
the photon frequency (the particle velocity) before and 
after the collision, 9 and (p are the scattering angles. 
There are five unknown parameters in the Compton 



scattering equations (34): ly, f3o, P,9,ip. If we fix (p 
and exploit (34), then we can obtain a single equation 
F{i', Pq) = that relates the frequency ly of the scattered 
light to the velocity c/3o of the particle. By linearizing 
such equation we obtain 

ly^A + BiPo-iM), (35) 

with A, B constants and c(/3o) mean value of the particle 
velocity. It is clear that, in order to fix the angle p, a 
coincidence system between the detector measuring the 
angle p of the scattered particle and the screen is needed. 

If the refractive index of the prism is a suitable and 
known function of the frequency v of the scattered light 
and assuming that linearization (35) holds, then we col- 
lect on a screen a spot whose location depends on the 
position-momentum correlation of the measured parti- 
cle. If without the prism the photon hits the screen at 
the position ^ = g, then with the prism ^ is displaced 
by an amount depending on the frequency v. If the re- 
fractive index is a linear function of ly, then we obtain 
^ = q + Cp, with C constant. Considering a more generic 
suitable nonlinear dependence of the refractive index on 
V, we obtain 

^^q + Cp + D{qp + pq), (36) 

with D constant. Since, as we have said above, (g) and 
(p) can be measured by means of other gedanken experi- 
ments, then the prism microscope is in principle suitable 
for the measurement of the correlator {qp + pq) . 

It is of course understood that experiments (i)-(iii) 
must be repeated many times, with identically prepared 
wave functions, to measure the above quantities (g), 
(Ag)^, (p), (Ap)^ and K with sufficient statistical ac- 
curacy. 

V. CONCLUSIONS 

We have shown that the Fermi gp function provides 
a very simple and nice phase space visualization of the 
dynamics of Gaussian wave packets. In particular, the 
area enclosed by the 5f = curve is of the order of the 
Planck constant and is a conserved quantity. We have 
also discussed a gedanken experiment for measuring the 
position-momentum correlation. This experiment adds 
to the well-known Heisenberg microscope and von Neum- 
man velocimeter, thus allowing the full determination of 
the state of a Gaussian wave-packet. 

The extension of the results obtained in this paper to 
more complex, non-Gaussian states has to face both nu- 
merical and conceptual difficulties. It is indeed clear from 
Eq. (5) that numerical errors in determining the gp — Q 
curve become relevant in the asymptotic region of large q 
where the wave function amplitude p is small. While for 
a Gaussian packet the gp = curve is sufficient to fully 
determine the state of the system, this is not the case 
for a generic and a priori unknown wave packet. The 
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complete determination of the generic state of a system 
requires consideration of the complex values of p±, ob- 
tained from Eq. (5) when dggp > 0. We then obtain 

= n^'-^ = -{^^)\ (37) 

from which the gp operator (3). and consequently the 
wave function ^(g, t) are determined. 

With regard to the comparison with well-know 
phase space distributions, notably the Wigncr function 
PwiQjPjty^'^^ , we first note that for Gaussian pack- 
ets the gp = curve coincide with the level curve of 
the points {q,p) such that pw{q,P,t) — ■^Piy^, with 
maximum value of pw- Different contour lev- 
els of pw correspond to different "equipotential curves" 
gp = constant. However, the gp function is conceptually 
very different from the Wigner function. First of all, not 
only gp fulfills Eq.(l) but also other operators such as 
gjr, 9%,— satisfy = 0, g'p'4) = 0,... . Ah informa- 
tion is enclosed in the zeroes of the Fermi function, from 
which one can reconstruct the wave function ijj. Finally, 
for the gp we cannot conceive any interpretation in terms 
of quasiprobabilities as for the Wigner function. 

The great simplicity of the Gaussian cases is due to 
the fact that the phase space region enclosed by the 
gp{q,p,t) = curve evolves in time as the correspond- 
ing classical phase space distribution pc{Q,P,t), uniform 
inside such region and zero-valued outside and whose 
dynamics is governed by the Liouville equation. In 



studying the free evolution of the superposition of two 
Gaussian wave packets we have seen that, as expected, 
this quantum-classical correspondence is broken. Due to 
quantum interference the gp = curve exhibits a rich, 
non-classical structure and, in particular, the area of the 
region enclosed by this curve does not remain constant 
over time. Therefore, quantum interference effects might 
be studied in phase space by means of the gp Fermi func- 
tion. 

Finally, we point out that it is possible to derive equa- 
tions governing the time evolution of the gp = curve^^. 
Such equations for generic, non-Gaussian wave packets 
differ from the classical Hamiltonian flow by the addition 
of fi-dependent terms. Therefore, one could use the Fermi 
function to study the quantum-classical correspondence. 
According to the Ehrenfest theorem, the propagation of 
a quantum mechanical wave packet is described for short 
times by classical equations of motion. The time at which 
this correspondence breaks down is called the Ehrenfest 
time^'^. After that time, due to quantum interference an 
initially Gaussian wave packet is "destroyed" , in that it 
splits into new small packets'*. Such splitting may be 
detected also by means of the Fermi function. 
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